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We present calculations for resonant inelastic x-ray scattering (RIXS) in edge-shared copper oxide 
systems, such as Li2Cu02 and CuGeOs, appropriate for hard x-ray scattering such as the copper 
if-edge. We perform exact diagonalizations of the multi-band Hubbard model and determine the 
energies, orbital character and resonance profiles of excitations which can be probed via RIXS. We 
find excellent agreement with recent results on I^CuCb and CuGeOs in the 2-7 eV photon energy 
loss range. 
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I. INTRODUCTION 

Techniques that probe the momentum dependence of 
strongly correlated electrons provide valuable informa- 
tion that can be used as a test for models appropriate 
for the cuprates. Angle-resolved photoemission (ARPES) 
provides important information about single-particle ex- 
citations of occupied states. 1 ' 2 Electronic Raman scatter- 
ing reveals information about multiparticle excitations 
at long wavelengths. 3 Resonant inelastic x-ray scattering 
(RIXS), a form of Raman spectroscopy, recently has un- 
dergone vast improvements, allowing probes of momen- 
tum dependent excitations into unoccupied states. 4 Im- 
provements in resolution now offer the possibility of a 
direct comparison of ARPES, Raman, and RIXS spectra 
in the several eV range with the hope of understanding 
the nature of strong correlations in the cuprates, and 
high temperature superconductivity, by studying elec- 
tronic excitations at several energy scales. 

As a fundamental building block of low-energy effective 
theories of the cuprates, the Zhang-Rice singlet (ZRS) 5 
is the starting point of the single-band Hubbard and t-J 
models widely believed to capture the low energy physics 
of strongly correlated systems. It is well-known that the 
ZRS features prominently in RIXS in the cuprates. A 
particular test of models of the cuprates lies in the differ- 
ence between RIXS in edge-shared versus corner-shared 
copper oxide systems. In materials such as La2Cu04 
and other corner sharing 2D members of the high T c 
family, the ZRS is stabilized by a strong gain in anti- 
ferromagnetic exchange energy between Cu and O, and 
can propagate to neighboring CUO4 plaquettes, with ef- 
fective hopping t = —0.35 eV along the Cu-0 bond di- 
rection, and t' = 0.15 eV along the diagonal. 5 In edge- 
shared cuprates, such as Li2Cu02 and CuGeOs, the Cu- 
O-Cu bond approximately forms a 90°-angle, 6 compared 
to 180° in corner-shared cuprates, and therefore accord- 
ing to Goodenough-Kanamori- Anderson, 7 the nearest- 
neighbor copper-copper exchange tends to be slightly fer- 
romagnetic. However, higher order processes like Cu-O- 
O-Cu hopping also contribute comparable antiferromag- 
netic exchange terms. For these reasons, small bond- 



angle variations among different edge-sharing compounds 
can make the exchange slightly anti-ferromagnetic or fer- 
romagnetic. In either case, the ZRS is less stable and 
mobile due to hybridization among different oxygen or- 
bitals not aligned with Cu d x i_ y i orbitals. 6 

Recent Cu K— edge RIXS data on edge-sharing 
cuprate Li2Cu02 8 ' 9 has revealed a small, non-dispersive 
peak at 2.1 eV, attributed to interatomic d—d excitations, 
a strong peak at 5.4 eV, and a weak peak at 7.6 eV, both 
attributed to charge transfer excitations. The peak at 5.4 
eV is doubly resonant for incident photon energies near 
8986 and 9000 eV. 8 In CuGeOs, clear excitations were 
found are 3.8 eV (weak) and 6 eV (strong), 9 ' 10 ' 11 and 
more recently 1.7 eV, 12 which were attributed in Ref. 12 
as a intra-atomic d — d excitation (1.7 eV), a ZRS ex- 
citation (3.8 eV), and a charge transfer excitation on a 
single copper oxide plaquette (6.4 eV). Similar features 
have also appeared in oxygen Is RIXS. 13 While a few 
momentum-dependent studies have been performed, the 
incident photon energy dependence, or resonance profile, 
is usually shown only for a few selected points in the 
Brillouin zone. 

In this paper we present a theory for RIXS for the cop- 
per if -edge and present exact diagonalizations of Hub- 
bard clusters of edge-sharing copper oxide plaquettes 
shown in Fig. 1. We calculate both the RIXS spectrum 
and the resonance profile for several cluster geometries, 
and determine the nature of excitations accessible via 
light scattering. We pay specific attention to the inter- 
mediate core-hole and final excited-state wavefunctions of 
RIXS accessible states. We find good agreement with the 
results on edge-shared copper oxide systems, and confirm 
a number of previous excitation assignments. 10 ' 11 ' 12 



II. MODEL AND METHOD 

A. RIXS process 

The RIXS process consists in scattering an incident 
photon of energy Hlji off a sample in the ground state 
|^o) having energy Eq. For if —edge RIXS, an electron 
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from the Is copper-core is photoexcited into the Cu 4p- 
band, leaving behind a core hole. The strong Coulomb 
repulsion between the core hole and d x 2_ y 2 holes on the 
same Cu atom reorganizes charge density, forming an in- 
termediate state \ipci) of eigenenergy E C {. The Ap electron 
resides in a rather extended wave-function and therefore 
interacts weakly with d and core states. We neglect its 
interaction hence providing us with a rather simple de- 
scription of RIXS. The Ap electrons then enter into the 
problem as a spectator at energy e± p = Huji — (E ci — Eq) 
and the intensity will be proportional to the projected 



4p-band density of states (DOS) at that energy. The in- 
cident polarization determines the necessary orbital pro- 
jection. The 4p electron then recombines with the Is 
core hole emitting a photon of energy hup, leaving the 
system in the final state \tjjf) of eigenenergy Ep and a 
photon energy transfer HQ = Huji — hwp. This leads to 
a simplification of the general expression given in Ref.4. 
In this approximation the spectrum will be given by the 
following relation and with the hujj dependence given by 
a convolution with the 4p band DOS : 
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Here C4 P -i s represents the Cu Is — 4p energy separation. 
For simplicity, in a first approach, we neglect specific Ap 
orientation and photon polarizations. The polarization 
dependence, as well as the related 4p orientation and the 
4p projected DOS, will be discussed in the last part of 
Section III. The parameter Ti represents damping of the 
intermediate state due to fluorescence and Auger Is de- 
cay, which we take as 1.0 eV. While this value determines 
the overall resonant enhancement, it does not affect the 
resolution of the energy-loss peaks in the spectrum. We 
broaden the delta function with a width O.leV to rep- 
resent current instrument resolution. We note that in 
Eq.(l) the states \^ C i) are identical to those in a core 
x-ray photoemission experiment. 

The Hamiltonian describing the open boundary cluster 
shown in Fig. (1), with rotated x and y local-directions 
along the p orbitals, is given by TL = Y^^^fja + 

£<(£<, «f +HY) + Uq E^, [1 - s$ J, 8 0f<r ,], with 

H-ija = tpdPij Pi yC r d j,(T ~\-tpd z Pij Pi j(T dz, j,a 
-\-t pp P i j Pi y(T Pj,<j + t pp p\ yCr Pj,a + h.C 

H\ = e p (pl a Pi,a + pl'aPi,*) + e dd[ a di^ + e dz d\. a d Zi a 
= Up P pl^PiApl^i + Up P pl^Pi^Pi d Pi,i 

+U dd d\ A d iA d\ d di^ + U d d d i iA d ziA d ii,i d zi,i- 

Here d\ a , d\ i (j ,p\ a ,p\ a and di^,d Zi a ,p iy(J ,p' i (J creates 
and annihilates a d x 2_ y 2, d 3z 2_ r 2, planar p or p' hole 
with spin a at site z, respectively, and the sum runs over 
nearest neighbors. The hopping amplitude t pp refers to 
the hopping between two nearest p x p y orbitals, and t 1 
refers to the hopping between two nearest-neighbor p a 
oxygen orbitals, where a = x, y. The overlap phase fac- 
tors are chosen such that Pij = 1 for neighbors in the 
x, — y direction, and -1 for — x, y. The oxygen phase fac- 
tors are analogously defined as P[ - = 1 for oxygens at 



x, y and — x, — y, and -1 for the other combinations. Here 
Sq (resp. so) creates (annihilates) a Cui s core-electron, 
and Uq is the core - d hole repulsion. For simplicity we 
neglect interactions involving the 4p photo-excited elec- 
tron, such as the 4p — 3d Coulomb interactions. These in- 
clude off-diagonal couplings between d states which may 
be important for singlet-triplet energy differences as well 
as providing pathways for symmetry- allowed d — d exci- 
tations. We also neglect exchange interactions involving 
the Is core-hole, which are very small. 

In what follows, we use a standard set of parameters 
[in eV] 5 : t pd = 1.1, t pp = 0.b : t f pp = -0.538 x t pp ,e d = 
0, 6 dz 1-7, A = e p - e d = 3.5, U pp 6, U dd 8.8, and 
Uq = 8. These parameters give a value for the copper 
magnetic exchange J = 0.14 eV for corner-shared copper 
oxides, and 2.8 meV for edge-shared. 



B. Discussion of the method and comparison with 
previous works 

We first perform full exact-diagonalization of the 
Hamiltonian without the core-hole, and then redo the 
calculation with the core-hole. Since the Hamiltonian 
conserves the total spin, we diagonalized in the S=l/2 
sector. From the eigenvalues and eigenvectors of the two 
steps, the RIXS spectrum for a given huji is computed 
from Eq.(l). 

As our calculations involve determining the eigenvec- 
tor of each state used to construct the matrix elements 
involved in RIXS, we know the character of each excita- 
tion exactly and do not need to make a priori assumptions 
about which peak corresponds to which fundamental ex- 
citation. This allows us to make direct correspondance 
between RIXS and other core-level spectroscopic tech- 
niques, without any further assumptions. 

This is quite different from the Lanczos method, which 
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FIG. 1: Top: Edge-shared Cu 5 0i 2 cluster. Bottom: XPS 
spectrum as a function of photon energy transfer and inci- 
dent photon energy for corner and edge-sharing copper-oxide 
clusters. 



truncates the Hilbert space and focuses on the low-lying 
energy excitations. Whereas this method allows investi- 
gation of slightly larger clusters, extensions to higher en- 
ergies is problematic particularly for RIXS where eigen- 
vectors are needed to evaluate the matrix elements be- 
tween higher energy states in order to obtain both the 
resonance profile (dependence on incident photon ener- 
gies) as well as the intensities of prominent peaks in the 
spectrum. Since the core-hole interaction is essentially 
in nature (on-site Coulomb repulsion), leading to only a 
slight derealization of the charge, we checked that the 
cluster-size does not play a crucial role in the qualitative 
shape of the spectra. In that sense, the calculations done 
in the present paper add more information and provide a 
more direct access to the relevant RIXS excitations com- 
pared to other techniques used in previous studies in the 
context of standard cuprates. 14 These studies treated the 
one-band Hubbard model where the Zhang-Rice singlets 
are, by construction, already formed, and focused only on 
the momentum dependence of the Mott gap excitation. 
By treating the problem in a multi-band model, the for- 
mation of Zhang- Rice excitations comes naturally in the 
problem and its signature in the RIXS intensity, peak dis- 
persion and resonance profile can be directly addressed. 
In fact it is widely recognized that the gap is a charge 
transfer gap which involves transition from oxygen 2p to 
copper 3d states. 15 



charge transfer (CT) excitation (d 10 L) at 3.5 eV whereby 
the screening of the core hole comes from the hole transfer 
from copper to oxygen on the same plaquette, 16 and 3) a 
10 eV high-energy poorly screened d 9 Cu hole on the core- 
hole site. These results are quantitatively similar to those 
presented earlier, 17 apart from a more realistic value of 
the exchange J in the corner-shared geometry. We note 
that x-ray absorption spectra in an extended system can 
be obtained from our XPS spectrum via a convolution 
with the copper Ap DOS, giving a broad absorption edge 
followed by spectral intensity extending higher in energy 
governed by the 4p bandwidth. 

These excitations, present in both corner and edge- 
shared geometries, occur at roughly the same energy 
since they are governed by charge transfer energies and 
the exchange oc t^ d / A between copper and oxygen: the 
Cu-Cu exchange does not play a major role in setting 
these energy scales. 

The main difference concerns the prominence of the 
non-locally screened d-intersite ZRS excitation in the 
corner-shared compared with the edge- shared systems. 
While the ZRS accounts for nearly 80 percent of the 
ground state wavefunction including the core hole for 
the corner-shared plaquettes, it is approximately 56 % 
for the edge-shared system, reflecting the number of oyx- 
gen hybridization pathways available in the edge-shared 
configuration. 18 Thus while the energy is comparable in 
both systems, the matrix elements are much weaker at 
the ZRS scale in the edge-shared compared to corner- 
shared systems. In addition, the mobility of the ZRS 
for each geometry determines the peak widths, reflecting 
the splitting between bonding and anti-bonding ZRS on 
sites neighboring the core hole. In the case of the corner- 
shared systems, this corresponds to 2t = — 0.7eF, while 
due to the orientation of the orbitals in the edge-shared 
systems, the splitting is 2t' = 0.3 eV. As a result, the 
locally screened d 10 L peak at 3.1 eV is much larger in 
edge-shared compared to corner-shared. This has led the 
phenomenon of high-temperature superconductivity, oc- 
curring in corner-shared compounds and absent in edge- 
shared compounds, to be associated with the stability of 
the ZRS. 19 We note that the primary role of the oxygen 
orientation in the formation of the ZRS is not captured 
in single-band Hubbard model calculations. 



III. RESULTS AND DISCUSSIONS 

A. Preliminary remarks 

In core level x-ray photoelectron spectroscopy (XPS), 
the final state corresponds to an intermediate state in 
the RIXS process. We first plot the XPS spectrum for 
both edge-shared cluster CU5O12 and corner-shared clus- 
ter CU5O16 in Fig. (1). In both cases, the prominent 
features correspond largely to 1) a non-locally screened 
d-intersite ZRS (d 9 L) feature at 2.2 eV, where the core 
hole has pushed a d 9 hole from the on-site Cu onto the 
oxygen ligands in the neighboring plaquette, 2) a local 



B. RIXS results : Attribution of the peaks 

The aim of this subsection is to discuss the overall 
shape of the RIXS spectrum and resonant profile for a 
CU3O8 cluster in the inset of Fig. 2a. We show how to 
attribute to each peak the corresponding excitation, as 
an example, we have chosen the same parameters as for 
XPS in Fig. 1, appropriate for Li2Cu02- From here on 
we set h equal to 1. Here the intermediate state core- hole 
is included on the edge of the cluster. The calculations 
are also performed for different locations of the core hole 
on the cluster, which has one hole per copper-plaquette. 
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FIG. 2: a) RIXS spectra for C113O8 cluster as a function 
of photon energy transfer Q, and incident photon energy uji, 
measured relative to the binding energy 64 p -i s . The inset 
shows a contour map over large frequencies. Matrix elements 
for CU3O8 are shown in b) & c), for the main peaks in RIXS at 
(ui - e 4p -i s ,n) = (2.2,2.2) and (3.1,5.07), eV respectively. 
Here the left, right eigenvector is the core hole, final state, 
respectively, and the dot and the circle represent a singlet 
state formed by a hole on the copper site (dot) and a hole 
delocalized on the oxygens having the ZR symmetry (circle). 



While the results differ slightly for different locations of 
the core hole, we note that the changes are only quanti- 
tative and give different broadening of the RIXS peaks 
due to the different mobility of the d-intersite Zhang- Rice 
excitation on the end or middle of the cluster. However, 
since the aim of this paper is to compare trends for dif- 
ferent cluster geometries and charge transfer energies, we 
will not discuss further these relatively minor changes in 
the spectra. 

The RIXS spectrum shown in the inset of Fig. 2a 
consists of two prominent excitations at approximately 
Q = 5.07 eV energy transfer, at incident photon energies 
of 3.1 and 10 eV. This is the local CT excitation, res- 
onant at both the locally screened d 10 L CT and poorly 
screened d 9 intermediate states, as shown in Fig. 1. In 
addition a weaker d-intersite ZRS-derived peak is present 
at fi = 2.2 eV, resonant for incident photon energies near 
2.2 eV above the Is — 4p transition. The XPS spectra 
may be crudely viewed as a sum over all Raman shifts of 
the RIXS spectrum as a function of the incident photon 
energy huoi. 

In Fig. 2a we show a blow up of the low energy RIXS 
spectrum. By examining the symmetries of the core-hole 
and final state wavefunctions, the nature of the RIXS ex- 
citations, and the origin of their resonant profile, can be 
determined by inspection of matrix elements appearing 
in Eq. 1. The ground state without the core hole | ^0) 
is largely of d 9 character in each plaquette. Only those 
intermediate states which have a large overlap with the 
ground state are relevant to the RIXS process. To deter- 



mine which final states | Vvo) are probed by RIXS and at 
what resonant energy they occur, one inspects the over- 
lap of the core hole states | ip c ) reachable via the ground 
state without the core hole | ^0) , with other excited 
states I ipi^o) m the same Hilbert space. For example, 
since the ground state of the cluster | ^0) is a singlet, it 
cannot be coupled to triplet intermediate states having 
a core hole. We immediately conclude that for our clus- 
ter all triplet states do not appear in the RIXS process 
without spin-orbit coupling. 

We sketch schematically in Fig. 2b and 2c the lead- 
ing order contributions to the wave functions of the in- 
termediate state with the core-hole and the final state 
with large matrix overlap, which correspond to the main 
excitations shown in the RIXS spectra. While the in- 
termediate state is predominantly one configuration, the 
final states can either be bonding or anti-bonding combi- 
nations of excitations on the plaquettes, at energies ap- 
proximately 2 eV above Eo, split by ZRS hopping. 

In Fig. 2b, it is clear that the resonance peak at 
uo — €4p-i s = 2.2eV and energy transfer Q = 2.2eV comes 
from the combination of intermediate and final states cor- 
responding to the creation of the ZRS in the final state, 
having strong overlap with the non-local, well-screened 
d 9 L intermediate state. The two final states of bond- 
ing and anti-bonding combinations of ZRS on the two 
plaquettes couple to the intermediate state, with a weak 
splitting in edge-shared systems. Even with our broad- 
ening of O.leV the separate peaks are difficult to resolve. 

The final states of the higher energy peak correspond 
to bonding and anti-bonding combinations of <i 10 L, hav- 
ing strong overlap with the intermediate core hole state 
with the ligand on the core hole site as shown in Fig. 2c. 
This lies at higher energy uji — e^ v -\ s = 3.07eF, and is 
the most prominent final state in XPS shown in Fig. 1. 

We note that we do not find any prominent d — d ex- 
citations lying in the relevant energy range below 2 eV, 
even when including other d orbit als. 20 This follows from 
a symmetry analysis of the states, whereby d 9 L and d 10 L 
configurations cannot hybridize other Cu d orbitals and 
combinations of oxygen ligands. Thus direct d — d exci- 
tations are inaccessible in our model, in agreement with 
electron energy loss cluster calculations. 21 It is plausible 
that d — d excitations may arise if symmetry-breaking 
interactions with the 4p are included. This is differ- 
ent for oxygen Is RIXS, as the creation of the core hole 
breaks this symmetry, allowing direct d — d excitations, 
as determined in recent cluster calculations. 22 Moreover, 
d — d excitations may arise from multiple couplings in 
non-resonant scattering 23 which are not treated in our 
calculation. 



C. Influence of the 4p-States — Polarization 
dependence — Material dependence 

So far, our analysis did not include the copper-4p 
states. However, in order to obtain results more closely 
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connected to experiment and to include a photon polar- 
ization dependence for the Cu K-edge process, the Cu-4p 
DOS must be included into the calculations. Indeed, the 
Cu-ls core-electron can only be excited into the Ap- levels 
if there is a finite DOS for unoccupied states. Due to 
crystal field effects, the 4p DOS for different £p x ,y,z or- 
bitals provide a polarization dependent shift of the reso- 
nance energies and secondary satellites governed by the 
peaks of each Ap projected DOS as well as the overall 
bandwidth. The orientation of the incident and emitted 
photon polarizations determine which Ap projected DOS 
is accessed. 



Since the 4p-states are rather extended orbitals their 
interaction with the Hubbard subbands is quite weak. 
However, depending on the incoming light-polarization, 
matrix elements are selected so that one can access differ- 
ent linear combinations of Ap Xl Ap y and Ap z states. Since 
the DOS distributions are differents, the RIXS spectra re- 
flect these characteristics. Thus to deduce the influence 
of the 4p-states, we neglect minor real-space variations of 
the combined cluster-4p wavefunction and compute the 
convolution of our raw-RIXS spectrum with the calcu- 
lated Ap projected DOS N^e^p): 
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The projected Ap DOS is calculated using linearized 
augmented plane wave method (LAPW) in a coordinate 
system consistent with experiment. 24 In the following 
subsection we present results for different edge-shared 
cuprates Li 2 Cu02 and CuGeOs. 



1. Li 2 Cu0 2 

For Li2Cu02, the Ap projected DOS is shown in 
Fig. 3(a). The Ap z DOS has two sharp peaks at roughly 
17eV, with subdominant peaks as shoulders to the main 
peak . The resulting RIXS spectrum for z photon polar- 
izations shown in Fig. 3(b) thus has a shifted d-intersite 
ZRS resonance at e^ p — uoi ~ 19eV instead of 2.2eV in 
Fig. 2, yet remains sharply peaked even though the Ap z 
bandwidth is over 20 eV wide. The d 10 L CT and poorly 
screened excitation at £1 ~ 5eV are likewise shifted to 
higher resonance frequencies. 

The Ap x DOS in Fig. 3(a) has a double-peak structure, 
with a single sharp peak at 4 eV and a tightly structured 
set of peaks around 10 eV. The resulting RIXS spec- 
trum for photon polarizations along the x-axis, given in 
Fig. 3(c), clearly shows two split resonances at ~ 6 — 7eV 
and ~ 13eV, a shift of the raw RIXS spectrum (Fig. 2) 
by the 4eV and 10eF Ap x DOS peaks. The difference in 
the relative intensities between Fig. 3(b) and Fig. 3(c) is 
due to the fact that the DOS for the p z orbitals is more 
intense at its maximum than for the p x orbital. 

Fig. 4 displays a side-by-side comparison of the present 
theory and recent experimental data by Kim et al. 8 on 




FIG. 3: (a) Cu-Ap partial density of states for Li2Cu02. (b) 
RIXS spectrum for parameters corresponding to Li2Cu02, 
convoluted with the p z -DOS. (c) RIXS spectrum for parame- 
ters corresponding to Li2Cu02, convoluted with the p x -DOS. 
For both (b) and (c) we took a core-hole lifetime of T = 1.5eV 



Li2Cu02- The theory and experimental curves show 
striking similarities in both the energy transfer and the 
incoming photon energy range. The energy scale, promi- 
nence of the local CT excitation, and the double reso- 
nance profile shown in Figs. 1,2,4 correspond well to the 
RIXS experiments. We remark that while Ref. 8 at- 
tributed the weak 2.1 eV peak to an interatomic d — d 
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FIG. 4: (a) Experimental data for L12C11O2 by Kim et al. in 
Ref.8. (b) Calculated RIXS spectrum for parameters corre- 
sponding to Li2Cu02 convoluted with the 4p-DOS. 
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FIG. 5: RIXS spectrum as a function of photon energy trans- 
fer and incident photon energy for parameters corresponding 
to CuGeOs, as defined in the text. 



excitation, and the strong peak at 5.4 eV peak to lo- 
cal CT excitation, in view of our results we would con- 
clude that the low energy excitation rather corresponds 
to the d-intersite ZRS excitation. The weak peak ob- 
served at 7.6 eV however has no correspondance in our 
calculations. This peak may be associated with a local 
CT excitation into states with different ligand symmetry 
orthogonal to the ZRS ligand states, which may become 
accessible once oxygen Ap — 2p interactions are included. 



2. CuGe0 3 

As pointed out in Ref.6, the large valency of Ge 4+ 
nearby the CuO planes increases the Madelung energy 
difference between the Cu and planar O sites in CuGeOs 
compared to LiCuC>2, leading to the values A = 4.9 and 
3.2 eV, respectively, for each system. Since this energy 
scale largely determines the energies of the d-intersite 
ZRS and the separation in energy to the local CT exci- 
tation, the RIXS spectrum and resonance profile should 
be largely different in the two cases. 

In a first approach, as we did for Li2Cu02, we present 
in Fig. 5 a non-convoluted RIXS spectrum for parame- 
ters corresponding to CuGeOs, which are the same set of 
parameters as before apart from the change t p d=l. 23 eV, 
and A = 4.9 eV. The larger value of A leads to several 
changes compared to the RIXS spectrum for L12CUO2 
parameters (see Fig. 2. a). First, we note that the d- 
intersite ZRS peak has moved higher in energy from the 
d 9 and core-hole ground states, lying at a photon energy 
transfer Q of 3.1 eV, and resonant at a higher incident 
photon energy uoi — e4 P _i s =3.1 eV. The bonding- anti- 
bonding splitting is smaller and the resonant peaks are 
more sharp. In addition, the CT (d 10 L) excitations move 
higher in energy, lying at uji— 64 P _i s = 3.9 eV and Q = 6.2 




b 



FIG. 6: CuGe0 3 structure. 



eV. This energy cost arises from the higher oxygen site 
energies as well as the reduction in exchange energy gain 
with the Cu spin. 

The CuGeOs structure in Fig. 6 is different from the 
one of Li2Cu02: the copper chains run now along the c- 
direction and they are tilted from the main- axis planes. 
We present in Fig. 7(a) the Ap projected DOS for 
CuGeOs. The coordinates system for which the calcu- 
lation has been done corresponds to the one presented in 
Fig. 6. As one can see, the result is strongly dependent 
on polarization: the density of Ap y is peaked around ~ 8 
eV, the 4^-density has its stronger weight around ~ 18 
eV and is much broader. 

The corresponding RIXS results are shown for pho- 
ton polarizations along p y and p z in Fig. 7(b) and (c), 
respectively. The overall shape of the spectra does not 
change much compared to the unconvoluted results (see 
Fig. 5) : one can still identify clearly a low-energy Zhang- 
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FIG. 7: (a) Cu-4p partial density of states in CuGeOs calcu- 
lated in WIEN2k. (b) RIXS spectrum for parameters corre- 
sponding to CuGeOs, convoluted with the p y -DOS. (c) RIXS 
spectrum for parameters corresponding to CuGeOs, convo- 
luted with the p z -DOS. 



FIG. 9: (Left) Experimental RIXS spectrum for CuGeOs with 
polarization dependence. (Right) Exact- diagonalizat ion re- 
sults obtained after convolution with the 4p projected DOS 
of Fig.7.(a) (see Subsection IIIC) 
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FIG. 8: (Left) Experimental RIXS spectrum for CuGeOs- 
(Right) Theoretical spectrum obtained for the same polariza- 
tion. 



Rice peak and a brighter d 10 L peak at higher energy as 
well as the poorly screened d 9 excitation peak. The main 
difference is a shift on the incident photon energy axis, 
of ~ 12 eV for p y and ~ 21 eV for p z . 

These results are compared to the recent experimental 
spectrum obtained by Hill et al 25 in Fig. 8 and 9. In this 
experiment a CT (d 10 L) peak at Q « 6.5 eV and the ZRS 
peak as well, closer to the elastic line, at about ~ 3.8 eV, 
were identified. These two observed peaks correspond to 
the two peaks seen in Fig. 5. 

In Fig. 8 we have used an energy-dependent damp- 
ing for the Lorentzian representing the delta function of 
Eq.(l). This energy dependence of the line width ac- 
counts for the fact that the higher energy states lie in 



a high density continua. For fr£t < 5.5eV, a quadratic 
energy dependence of the damping rate was taken and 
assumed to saturate at higher Raman shifts. This form 
was chosen in order to account for the relative intensity 
and width of the ZR and d 10 L peaks. 

We note that the slight discrepancy between theory 
and experiment in the position and intensity of the d- 
intersite ZRS may be remedied by fine-tuning the clus- 
ter parameters, to more accurately reflect the material 
properties (such as J and the effective ZR- hopping t), as 
well as incorporating a more realistic resolution broad- 
ening parameter. Another difference between theory and 
experiment is the relative intensity of the peaks for dif- 
ferent polarizations: while the peaks are more intense 
for a polarization along b in the experiment, it is the 
opposite for the simulation. This issue may be related 
to the two following points: (i) from a theoretical point 
of view the relative intensity of the peaks is sensitively 
linked to the relative peaks height of the DOS (since in 
Eq.(l) it appears to be squared), thus a slight variation 
in the DFT calculation might make a difference, (ii) Due 
to experimental uncertainties a direct comparison of the 
measured absolute intensities for e||b and e||c is not pos- 
sible. However, for a fixed polarization the measured 
relative intensities of the ZR and the CT feature can be 
compared. 

Finally we remark that the crystal structures of 
CuGeOs and Li2Cu02 are strongly different (see for in- 
stance Refs. 26,27). In both cases the Cu02 chains run 
along a principal axis of the crystal but in CuGeOs 
they are not oriented along the same plane, whereas for 
Li2Cu02 they are. This has strong implications con- 
cerning the interpretation of the experimental results : 
for CuGeOs the polarization-dependent experiments will 
give RIXS spectra which are a result of convolution with 
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the projected DOS in a rotated basis. The result of such 
a convolution is given in the right panel of Fig. 9. 

IV. CONCLUSIONS 

A full exact-diagonalization treatment of multi-band 
Hubbard clusters has been used to calculate the polar- 
ization, incident photon energy, and transfered photon 
energy dependencies, of RIXS in one-dimensional, edge- 
shared insulating cuprates. By building upon previous 
work, 16 the RIXS spectrum can be constructed using 
XPS final-states as the intermediate states of inelastic 
x-ray scattering. The intensity and resonance profile, as 
well as the identity and character of the main RIXS peaks 
has been obtained using the eigenstates and eigenenergies 
of clusters with and without the core-hole. These peaks 
correspond to non- locally screened d-intersite ZRS d 9 L, 
local CT d 10 L, and poorly screened d 9 , each having in- 
tensities and resonance profile determined by the overlap 
of these states with the intermediate core- hole states. 

The polarization dependent RIXS spectra were ob- 
tained by convolving the raw RIXS spectra with the cop- 
per Ap projected DOS determined by WIEN2k density 
functional calculations. While the resulting RIXS spec- 
tra includes the Ap states as merely spectators to the 
re- arrangement of copper-oxygen valence charge density, 
a polarization and material dependent RIXS emerges due 
to the peak structures and intensities of the Ap orbitally- 
projected intermediate states. 

Direct comparison (see Figs. 4,9) with recent 
data for two edge-shared cuprates, Li2Cu02 and 
CuGeOs, 8,9,10,11,12 shows both qualitative and, more im- 
portantly, quantitative agreements between theory and 
experiment for the shape, energy, width, and resonance 
profile of RIXS spectral features. This indicates that the 
local cluster model can adequately reproduce the data 
for photon energy losses in the energy range of 2-7 eV. 
We again remark that d — d excitations are not cap- 
tured in the model, and further, that spin-flip excita- 
tions, which require large clusters to observe bi-magnon 
spin rearrangements at long wavelengths, are also very 
weak. Better agreement of the calculated spectra with 



experiment may be obtained by using more refined clus- 
ter parameters, multi-pole couplings, and more accurate 
representation of resolution broadening. 

A polarization dependent RIXS experimental study, 
for Li2Cu02, should find significant differences in the res- 
onant profiles for in-plane polarization and polarization 
perpendicular to the Cu02-plane. Fig. 3 indicates that 
peak intensities are brighter and their resonance ener- 
gies are shifted by ~ 10 eV for in-plane compared to 
perpendicular polarizations. These predictions as well as 
others, like multiple resonant features resulting from 4p- 
projected states, open a way to infer from Cu K— edge 
and other indirect RIXS experiments the unoccupied 
DOS. In this way RIXS can be seen as a complemen- 
tary tool to core-level spectroscopies such as XPS and 
XAS, providing at the same time diverse and impor- 
tant spectral information on elementary excitations in 
strongly correlated materials. 

In summary, we have explored the Cu K— edge spec- 
trum for different edge-shared copper oxide systems, and 
have shown that the energy scales, peak intensities, and 
resonance profiles can be explained qualitatively within 
the context of exact diagonalization studies of the multi- 
band Hubbard model. The prospect of mapping out the 
full momentum, and resonance profile dependences would 
allow for a direct spectral check of the many roles played 
by local electronic correlations in the multi-band Hub- 
bard approach. This remains a topic of future research. 
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